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Abstract 

Microcanonical description is characterized by the presence of an internal symmetry closely related 
with the dynamical origin of this ensemble: the reparametrization invariance. Such symmetry pos- 
sibilities the development of a non Riemannian geometric formulation within the microcanonical 
description of an isolated system, which leads to an unexpected generalization of the Gibbs canon- 
ical ensemble and the classical fluctuation theory for the open systems, the improvement of Monte 
Carlo simulations based on the canonical ensemble, as well as a reconsideration of any classification 
scheme of the phase transitions based on the concavity of the microcanonical entropy. 
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Generally speaking, a phase transition brings about a 
sudden change of the macroscopic properties of a sys- 
tem while smoothly varying of a control parameter. The 
mathematical description of phase transitions in the con- 
ventional Thermodynamics is based on the loss of analyt- 
icity (the appearance of a singularity in the firsts deriva- 
tives) of the (gran) canonical thermodynamic potential 
(0; 0), which are related with the existence of zeros in 
the partition function appearing with the imposition of 
the thermodynamic limit (|3|). 

However, phase transitions are also relevant in systems 
outside the thermodynamic limit. A nontrivial example 
is the nuclear multifragmentation, a phenomenon taking 
place as a consequence of peripheral collisions of heavy 
ions, which can be classified under the light of the new 
developments as a first-order phase transition (0;0). 

The key for taking into account the phase transitions in 
isolated systems outside the thermodynamic limit is the 
consideration of the well-known Boltzmann Principle: 



]xlW 



(1) 



his celebrated gravestone epitaph. Since W is the num- 
ber of microscopic states compatible with a given macro- 
scopic state, the Boltzmann entropy 2] is a measure of 
the size of the microcanonical ensemble. Within Clas- 
sical Statistical Mechanics W is just the microcanonical 
accessible phase space volume, that is, a geometric quan- 
tity. This explains why the entropy ^ does not satisfy 
the concavity and the extensivity properties, neither de- 
mands the imposition of the thermodynamic limit or a 
probabilistic interpretation like the Shannon-Boltzmann 
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Gibbs extensive entropy: 



k 



(2) 



As already shown by Gross (0; 0), phase transitions in 
small systems could be classified within the microcanon- 
ical ensemble throughout the "topology" of the Hessian 
of the Boltzmann entropy. 

We will show in the present work that the microcanon- 
ical description is characterized by the existence of an in- 
ternal symmetry: the reparametrization invariance. We 
shall demonstrate that the presence of this symmetry im- 
plies a revision of classification of phase transitions based 
on the concavity of the entropy and the theory of the sta- 
tistical ensembles derived from microcanonical basis. 



II. REPARAMETRIZATION INVARIANCE 
A. Geometrical basis 

Universality of the microscopic mechanisms of chaotic- 
ity provides a general background for justifying the nec- 
essary ergodicity which supports a thermostatistical de- 
scription with microcanonical basis for all those noninte- 
grable many-body Hamiltonian systems Ipl Ifj fTTll: Hit Il2l . 
Thus, the microcanonical ensemble: 

is just a dynamical ensemble where every macroscopic 
characterization has a direct mechanical interpretation. 
Here, X represents a given point of the phase space T 

and I (X) = [i 1 (X) ,I 2 (X),... I" {X)\ are all those 

relevant (analytical) integrals of motion determining the 
microcanonical description in a given application (since 
Poincare- Fermi theorem Q0): the total energy, the an- 
gular and linear momentum). 

The admissible values of the set of integrals of motion 
I (X) could be considered as the "coordinate points" / = 
{l 1 , I 2 , . . . J™} of certain subset Hi of the n-dimensional 
Euclidean space H n . Each of these points determines 
certain sub-manifold S p of the phase space T : 



x e s p = {x e r |vfc i k (x) = i k } 



(4) 



in which the system trajectories spread uniformly in ac- 
cordance with the ergodic character of the microscopic 
dynamics. Such sub-manifolds defines a partition 3 of 
the phase space T in disjoint sub-manifolds: 



c r 



\Js p = r ; s p n s q = I 

v ) 



(5) 



Definitions 0] and [5] allow the existence of a bijective map 
ipi between the elements of 9 (sub-manifolds S p C T) 
and the elements of Hi (points / £ H n ): 
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Hi = {VS P e 3 (r) ]/eK ; c H n } . (6) 



Thus, the partition 5 has the same topological features 
of the n-dimensional Euclidean subset Hi. For this rea- 
son 3 will be referred as the abstract space of the inte- 
grals of motions. We say that the map ipi defines the 
n-dimensional Euclidean coordinate representation Hi of 
the abstract space 3. 

Let us now to consider another subset H v C H n with 
the same diffeomorphic structure of the subset Hi and 
the following diffeormorphic map tp among them: 



tp : Hi — > H^ 



VI £TZi 3(p£ H, 



det 



^0 



(7) 

We say that the map tp represents a general 
reparametrization change of the microcanonical descrip- 
tion since it allows us to introduce another n-dimensional 
Euclidean coordinate representation Hip by considering 
the bijective map ip v = ipiotp^ 1 : 



H u 



{VS P e 9= 3<p e H v c H n } . (8) 



The above reparametrization change tp also induces the 
following reparametrization of the relevant integrals of 
motion ipx '■ I {X) — » tp (X), where: 

<p (x) = [p 1 (i (x)) ,p 2 (/ (x)) , . . . v n (i (X))}. 

(9) 

Since / (X) are integrals of motions, every <p k (X)^ £ 

tp (X) will be also an integral of motion. The bijective 
character of the reparametrization change tp : Hi — ► H v 

allows us to say that the sets tp (X) and I (X) are equiva- 
lent representations of the relevant integrals of motion of 
the microcanonical description because of they generate 
the same phase space partition Q\ 

The interesting question is that the microcanonical 
ensemble is invariant under every reparametrization 
change. Considering the identity: 



S(tp-tp(X)) 



dip 
~dl 



S(I-I(X) 



(10) 



where \dtp/dl\ ^ is the Jacobian of the reparametriza- 
tion change tp, the phase space integration leads to the 
following transformation rule for the microcanonical par- 
tition function: 



n(tp,N) 



dp 
Hi 



n(i,N), 



(ii) 



leading in this way to the reparametrization invariance 
of the microcanonical distribution function: 



1 



n(tp,N) 



5(tp-tp(X)) = 



1 



ft (I, N) 



S(I-I(X)). (12) 



A corollary of the identity El is that the Physics 
derived from the microcanonical description is 
reparametrization invariant since the expectation 
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values (O) of any macroscopic observable O (X) ob- 
tained from the microcanonical distribution function 
u>m (X) exhibits this kind of symmetry: 

(O) = I 6 (X) lum (X) dX (O) (ip, N) = (O) {I, N) . 



(13) 

The reparametrization invariance does not introduce 
anything new in the macroscopic description of a given 
system, except the possibility of describing the micro- 
canonical macroscopic state by using any coordinate rep- 
resentation of the abstract space 9, a situation analogue 
to the possibility of describing the physical space 1Z 3 by 
using a Cartesian coordinates (x, y, z) or a spherical co- 
ordinates (r,9,tp). Thus, we can develop a geometrical 
formulation of Thermostatistics within the microcanoni- 
cal ensemble. 



B. Covariant transformation rules 

The microcanonical partition function allows us to in- 
troduce an invariant measure d/i = fldl for the ab- 
stract space 5, leading in this way to an invariant def- 
inition of the Boltzmann entropy Sb = l~nW, where 
W = J s dfi characterizes certain coarse grained par- 
tition {Sq, I U Q £a = 3}- In the thermodynamic limit 
TV — > oo the coarsed grained nature of the Boltzmann 
entropy can be disregarded and taken as a scalar func- 
tion defined on the space (see the appendix [BJ) 

Since the Boltzmann entropy is a scalar function, its 
first derivatives obey the transformation rule of a covari- 
ant vector during the reparametrization changes: 



dS B dl s dS E 



dip K df K dl s 



(14) 



However, the second derivatives (Hessian) of the entropy 
do not correspond to a second rank covariant tensor: 



d 2 S B 



di s di T d 2 s f 



d 2 I s dS f 



d<p K d<p L dtp K dip L 8I K dI T d<p K dtp L dl s 



(15) 



because of the correct transformation rule should be 
given by: 



dl s dI T 

d(p K dip- 



-TST- 



(16) 



Covariant tensors can be derived from the differentia- 
tion of vectors within a Riemannian geometry throughout 
the introduction of the covariant differentiation, which 
depends on the existence of an appropriate metric. There 
are in the past other geometric formulations of the ther- 
modynamics which identify the metric with Hessian of 
entropy. We provide in the appendix^ a little explana- 
tion of the different underlying physics supporting such 
geometrical constructions. A complete review of such 
formulations can be seen in ref. l)l3|) . 



The using of the Hessian as a metric will be discarded 
in the present approach due to the Hessian of the en- 
tropy associated to an isolated Hamiltonian system does 
not correspond to a second rank covariant vector. Reader 
may think that the non reparametrization invariance of 
the entropy Hessian is a ugly defect of the present geo- 
metric formulation. Contrary, we will show that this fea- 
ture could become a big advantage in understanding the 
nature of the phase transitions within the microcanoni- 
cal ensemble, as well as in the improving of some Monte 
Carlo methods based on the Statistical Mechanics. 

With some important exceptions like the astrophysi- 
cal systems, the nuclear, atomic and molecular clusters, 
and some other systems, most applications of statisti- 
cal mechanics are concerned with the behavior of atoms 
and molecules subjected to short-range forces of electro- 
static origin, i.e. gases, liquids, and solids. These sys- 
tems are typically enclosed by rigid boundaries, and con- 
sequently, the only relevant integral of motion here is the 
total energy E. It is very easy to verify that any bijec- 
tive application of the total energy Q = (E) ensures 
the reparametrization invariance of the microcanonical 
description: 
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n(6,JV) 



s(e-e 



N 



1 



ft (E, N) 



5(E-H 



(17) 



being 6^ = © (Hn^, and consequently, the results ob- 
tained above arc still applicable in this context. 

Summarizing the result of the present section: (1) The 
microcanonical description is reparametrization invari- 
ant; (2) The Boltzmann entropy is just a scalar function 
under the reparametrizations; (3) While the first deriva- 
tives of the entropy are the components of a covariant 
vector, the second derivatives do not correspond to a sec- 
ond rank covariant tensor. The Boltzmann entropy also 
allows the introduction of other thermodynamic relations 
which exhibit explicitly the reparametrization covariancc 
(see in appendix [UJ. 



III. THERMODYNAMIC AL IMPLICATIONS 

The present proposal comes from by arising the 
reparametrization invariance to a fundamental status 
within the microcanonical description. In the sake of 
simplicity, let us analyze the implications of this symme- 
try in the thermodynamical description of a Hamiltonian 
system with a microscopic dynamics driven by short- 
range forces whose microcanonical description is deter- 
mined only from the consideration of the total energy 
E, that is, an ordinary extensive system. 

The preliminary interest of the present analysis is how 
identify the phase transitions within the microcanonical 
description. Taking into account the dynamical origin of 
this ensemble, a phase transition within the microcanon- 
ical ensemble should be the macroscopic manifestation 
of certain sudden change in the microscopic level which 
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manifests itself as a mathematical anomaly of the Boltz- 
mann entropy. Since the entropy is a continuous scalar 
function, the most important mathematical anomalies of 
the entropy per particle s = Sb/N are the following: 

(A) Regions where s is not locally concave; 

(B) Every lost of analyticity in the thermodynamic 
limit N — ► oo. 

Anomaly type A is directly related with the first-order 
phase transitions in the conventional Thermodynamics, 
while we shall show that the anomaly type B can be as- 
sociated to the second-order phase transitions and other 
anomalous behaviors. We will concentrate in the next 
subsections to the analysis of the anomalies type A. 



A. Why is so important the concavity of the mi- 
crocanonical entropy? 

Most of applications of the Equilibrium Statistical Me- 
chanics start from the consideration of the Gibbs canon- 
ical ensemble: 



u c (0,N) = 



1 



exp 



{-p&n}, 



(18) 



which provides the macroscopic characterization of a 
Hamiltonian system with a thermal contact with a heat 
bath. This description becomes equivalent to the one 
carried out by using the microcanonical ensemble almost 
everywhere when N is large enough, which can be easily 
verified by considering the Laplace transformation be- 
tween their partition functions: 



Z (0, N) = / exp (-0E) n (E, N) dE. 



(19) 



This equation is conveniently rewritten by introducing 
the coarsed grained entropy S (E, N) — In {SI (E, N) Seq} 
and the Planck thermodynamic potential P (0, N) — 
-lnZ(0, N) 1 : 



Se'n 



-{/3E-S(E,N)} dE _ 



(20) 



The energy per particle dependence of the distribution 
function r\ (e; 0, N) of the interest system within the 
canonical ensemble is just the exponential function of the 
above integral: rj(e;0,N) cx exp (— N (0e — s (e, N))), 
being s(e;N) — S (Ne 7 N) /N the entropy per parti- 
cle. When N is large enough, r\ (e; 0) exhibits very sharp 



1 P ((3, N) is related with the Helmholtz free energy potential by 
F(]3,N) =/3- 1 P((3,N). 



peaks around all those stationary points {s m } minimiz- 
ing the functional p (e; 0, N) — (3e ~ s (e; N): 

/3 =^^and«=^%i^<0. (21) 



de 



de 2 



However, in spite of the existence of several peaks, 
the presence of the system size N in the argument of 
the exponential dependence of the distribution function 
rj (e; 0, N) allows the Planck potential per particle to take 
asymptotically the value of the global minimum e gm of 
p(e;(3,N): 



P (0 N) 

p (0, N) = ^ ' ~ 0e gm - s (e gm ;N) + O 



InN 



(22) 

which determines in practice the thermodynamic values 
of the microscopic observables within the canonical de- 
scription. 

There exist ensemble equivalence when there is only 
one peak at a given in the distribution function 
r\ (e; 0, N) (the macroscopic state at a given in the 
canonical ensemble is equivalent to the macroscopic state 
of the system within the microcanonical ensemble with 
£ = £gm)- The square dispersion of the energy (Ae 2 ) ~ 
— 1/ (Nk) here decreases with the increasing of the sys- 
tem size N, so that, the system energy is effectively fixed 
in the thermodynamic limit. When TV is large but finite, 
the equivalence is only asymptotic, that is, the canonical 
average (^^j of any microscopic observable A converge 
asymptotically towards the corresponding microcanoni- 
cal average (A) with the N increasing: 



6a 



A) - A 



1 

N' 



(23) 



The reader may notice that the equation 1221 is just the 
Legendre transformation P (/3, N) — (3E—S (E, N) which 
relates the thermodynamic potentials of the correspond- 
ing ensembles, whose validity justifies the applicability of 
the well-known thermodynamic formalism of the conven- 
tional Thermodynamics. 

On the other hand, there exist ensemble inequivalence 
when the distribution function rj (e; 0, N) exhibits a mul- 
timodal character, being this feature a direct manifesta- 
tion of the existence of several metastable "microcanon- 
ical states" within the canonical description at a given 
0. The energy interchange with the heat bath (thermo- 
stat) provokes here a random transition of the system 
energy among these metastable states, and consequently, 
the system undergoes the incidence of very large fluctu- 
ations and non homogeneities. 

Undoubtedly, this is the well-known phenomenon of 
phase coexistence usually referred as the occurrence of a 
first-order phase transition: any small variation of the 
thermostat parameter close to a critical value 0c pro- 
vokes that certain metastable state becomes in a global 
minima of the functional p (e; 0, N), leading in this way 
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to the occurrence of an abrupt change of the average en- 
ergy per particle (e) c of the system. Since the canonical 
average of the energy per particle can be derived from 
the Planck potential per particle p (/?; N) as follows: 



1 



dZ (0, N) _ dp {0; N) 



NZ (/?, N) 80 



80 



(24) 



such thermodynamic function exhibits a sudden change 
in its first derivative at the critical point 0c, which be- 
comes a discontinuity in the thermodynamic limit, being 
this behavior a feature of the first-order phase transitions 
in the conventional Thermodynamics. 

Since the transition from a given metastable state to 
another demands the incidence of very large fluctuations 
of the total energy, such events involve a very large char- 
acteristic time due to the probability of occurrence of a 
very large fluctuation decreases exponentially with the 
increasing of the system size N. This dynamical behav- 
ior leads to an exponential divergence of the correlations 
times with N during the Monte Carlo simulations based 
on the Gibbs canonical ensemble, a phenomenon referred 
as a supercritical slowing down |)18T) . 

Supposing the existence and continuity of the first and 
the second derivatives of the entropy, the existence of 
two or more stationary points satisfying the conditions 
EH! at a given is also associated to the presence of con- 
vex regions of the microcanonical entropy. Such convex 
regions represent thermodynamic states with a negative 
heat capacity within the microcanonical description: 



c (m) = f <^£_ 

dE 



- (ft*L\ fd 2 SY 
\8EJ \8E 2 ) 



(25) 



which is an anomalous behavior since the heat capacity 
in the canonical ensemble is always positive: 



C^ = ^ = 2 ((AEf 



> 0. 



(26) 



According to the conditions 1211 such anomalous ther- 
modynamic states will not be observed in the Gibbs 
canonical ensemble 1181 Particularly, the corresponding 
microcanonical entropy can not be obtained from the 
Legendre transformation S = 0E — P starting from 
the Planck potential P. Gross refers this situation as a 
cathastrophe of the Legendre transformation 10). Thus, 
the existence of an ensemble inequivalence associated to 
the convexity of the microcanonical entropy leads to a 
significant lost of information when the canonical descrip- 
tion is performed instead the microcanonical description: 
the Gibbs canonical ensemble is unable to describe the 
system features during the phase coexistence, i.e. the 
existence of a non vanishing interphase tension (0). 

Summarizing: The concavity of the Boltzmann en- 
tropy is a necessary condition for the ensemble equiv- 
alence and the validity of the Legendre transformation 
P = 0E — S between the thermodynamic potentials 



which supports the applicability of the well-known ther- 
modynamic formalism in the conventional Thermody- 
namics. The existence of convex regions in the Boltz- 
mann entropy of short-range interacting systems can be 
related with the existence of the first-order phase tran- 
sitions, which are also associated with the lost of ana- 
lyticity of the Planck P (0, N) or the Helmholtz free en- 
ergy F (0, N) potential (discontinuity of first derivatives) 
in the thermodynamic limit. The existence of the mul- 
timodal character of the canonical energy distribution 
function under the presence of first-order phase transi- 
tions is the origin of the phenomenon of critical slowing 
down in the neighborhood of the critical point during the 
Monte Carlo simulations based on the Gibbs canonical 
ensemble 1181 



B. The concave or convex character of the micro- 
canonical entropy is an ambiguous concept from the 
reparametrization invariance viewpoint! 

As already discussed in the subsection above, the con- 
cavity of the microcanonical entropy is a very impor- 
tant condition within the conventional Thermodynamics. 
However, the consideration of the reparametrization in- 
variance leads naturally to the ambiguous character of the 
concavity of the microcanonical entropy: the concave or 
convex character of a scalar function like the entropy de- 
pends crucially on the parametrization used for describe 
it. Let us consider a trivial example. 

Let s be a positive real map defined on a seminfinite 
Euclidean line £, s : C — > which is given by the con- 
cave function s (x) = \fx in the coordinate representation 
1Z X of £ (where x > 0). Let ip be a reparametrization 
change ip : 1Z X — ► 1Z V given by y — ip (x) — i* . The 
map s in the new representation !Z y of the seminfinite 
Euclidean line C is now given by the function s (y) = y 2 
(with y > 0), which is clearly a convex function. 

The ambiguous character of the microcanonical 
entropy is straightforwardly followed from the non 
reparametrization invariance of the entropy Hessian 
shown in the eauation ll5l The only exceptions are those 
thermodynamic states where the corresponding entropy 
shows a local extreme: the first derivatives of the en- 
tropy vanish there and the Hessian behaves eventually 
as a second rank tensor. An analogue situation was used 
by Ruppeiner for establish a Riemannian interpretation 
of the Thermodynamics in the ref. l(l^) (see in appendix 

Reader may object correctly that such ambiguity of the 
concave character of the entropy is irrelevant within the 
conventional Thermodynamics because of the only one 
admissible representation (E) in this framework is the 
total energy, 9 (E) = E. This viewpoint follows directly 
from the fact that most of applications of the Equilibrium 
Statistical Mechanics start from the consideration of the 
Gibbs canonical ensemble 1181 for more general, from the 
Boltzmann- Gibbs distributions) . However, it is necessary 
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to recall that the reparametrization invariance is only rel- 
evant in the microcanonical description since the Gibbs 
canonical ensembleEDis not reparametrization invariant. 
Their different background physical conditions are cru- 
cial: the microcanonical ensemble describes an isolated 
system while the Gibbs canonical ensemble describes the 
same system with a weak interaction with a thermostat, 
that is, an open system. 

At first glance, the above ambiguity suggests us that 
the first-order phase transitions associated with the ex- 
istence of convex regions of the entropy could not be 
a phenomenon microcanonically relevant due to the 
reparametrization invariance of this ensemble. This pre- 
liminary idea contradicts our common sense educated 
in the everyday practice where the phase coexistence 
is a very familiar phenomenon. Nevertheless, most of 
these experiences are observed under those background 
conditions which lead to the applicability of the Gibbs 
canonical ensemble El We shall demonstrate in the next 
subsection that this ensemble can be generalized in or- 
der to develop a thermodynamic framework where the 
reparametrization invariance plays a more fundamental 
role than in the conventional Thermodynamics. 

C. Generalized canonical ensemble 

The Gibbs canonical ensemble EJ belongs to a fam- 
ily of equilibrium distribution functions called the 
Boltzmann- Gibbs distributions (jlj). Generally speaking, 
the Boltzmann-Gibbs distributions provide the macro- 
scopic characterization of an open system in thermody- 
namic equilibrium with its very large surrounding, where 
this system and its surrounding constitute a closed envi- 
ronment. Of course, the specific form of this distribution 
functions depends on the external conditions which are 
been imposed to the interest system, i.e.: the well-known 
Gran canonical ensemble: 

fiOJ,M) = ^«p{^(A + ^)}, (27) 

describes a system which is able to interchange energy 
and particles with its surrounding (reservoir). 

Although the Boltzmann-Gibbs distribution functions 
play a fundamental role in most of applications of the 
Statistical Mechanics and the classical fluctuation the- 
ory, they are not the only physically admissible: an open 
system could be arbitrarily affected by many other exter- 
nal conditions so that the whole experimental setup (the 
interest system + its surrounding) do not correspond nec- 
essarily to a closed environment. It is important to men- 
tion that although an open system should exhibit a sta- 
tionary macroscopic state when the external conditions 
are also stationary, most of these stationary states do not 
correspond to equilibrium conditions since the probabil- 
ity currents between microstates may not vanish. 

Let us consider certain experimental setup which al- 
lows the interest short-range Hamiltonian system Hn to 



interchange energy with its surroundings in some station- 
ary way that leads this system to a thermodynamic equi- 
librium described by the following generalized canonical 
ensemble: 

Cj gc (v, N) = - 1 — exp (-A) , (28) 

where On = 9 (Hn^J corresponds to certain 

reparametrization of the total energy E — > 9 = 9 (E) = 
Nip (E/N) preserving the extensive character of the en- 
ergy in the thermodynamic limit. Obviously, the generic 
form of the distribution function contains the Gibbs 
canonical ensemble 1181 as a particular case, where the 
only explicit difference is the introduction of the energy 
reparametrization (E). 

The generalized distribution function |2U can be 
straightforward derived from the maximization of the 
Shannon-Boltzmann-Gibbs extensive entropy [21 by pre- 
serving the average value of the microscopic observable 
e(E): 

(0) =Y,®(Ek)Pk, (29) 

fe 

This observation allows us to understand that the role 
of the experimental setup which leads to the generalized 
canonical ensemble 1281 is just to preserve the average l2*9l 
instead the energy average: 

(E)=Y,EkPk- (30) 

fe 

Obviously, the present experimental setup is a more so- 
phisticate version of the thermostat of the Gibbs canoni- 
cal ensemble which keeps fixed the " canonical" parameter 
rj instead the inverse temperature (3. We shall refer this 
external setup as a generalized thermostat whose canon- 
ical parameter rj exhibits in the present framework the 
same thermodynamic relevance of the ordinary tempera- 
ture in the conventional Thermodynamics. 

The generalized canonical ensemble defines a complete 
family of equilibrium distribution functions which exhibit 
the same features of the Gibbs canonical ensemble al- 
most everywhere. This affirmation follows directly from 
the reparametrization invariance of the microcanonical 
ensemble: 

Q(E,N)dE = Sl(Q,N)d&, (31) 

which allows us to express the partition function of the 
generalized canonical ensemble (|28|) by using the follow- 
ing Laplace Transformation: 

Z (r), N) = J exp (-779) n (9, N) de. (32) 

The analogy between the equations |221 and ^| allows us 
to understand that the generalized canonical ensemble 
admits a simple extension of the features of the Gibbs 
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canonical ensemble by using the reparametrization O (E) 
instead the total energy E of the extensive system. 

Particularly, the generalized Planck potential 
P(rj,N) = —lnZ(r],N) in the thermodynamic limit 
iV — ► oo can be approximated by the global minimum 
®gm G {9 m } of the Legendre transformation: 



P(r),N)~ rje gm -S(@ gm ,N), 



(33) 



whose the stationary points {0 m } are derived from the 
conditions: 



dS(e m ,N) 9 2 S(0 m ,7V) 



90 



de 2 



< 0. 



(34) 



Thus, such generalized canonical description at a given r\ 
becomes equivalent in the thermodynamic limit N — > oo 
to a macroscopic state of the microcanonical ensemble 
when the global minimum Q gm is the only one station- 
ary point. We can refer this situation as a local ensemble 
equivalence between a given generalized canonical ensem- 
ble and the microcanonical ensemble. We can also refer 
to a global ensemble equivalence when these ensembles 
are locally equivalent everywhere. 

Two macroscopic states corresponding to different gen- 
eralized canonical descriptions (with reparametrizations 
0i (E) and ©2 (E) respectively) are mutually equivalent 
when these macroscopic states are locally equivalent to 
the same macroscopic state in the microcanonical ensem- 
ble. It is very easy to derive from the stationary condi- 
tions [Ml that the canonical parameters r)® 1 and n@ 2 of 
mutually equivalent macroscopic states are related by the 
following transformation rule: 



Ve 2 



dS 
90 



dS 
90i 



m 2 = 



901 
90' 



(35) 



This is exactly the transformation rule during the 
reparametrization changes of unidimcnsional covari- 
ant vectors expressed in the equation While 
a microcanonical state remains invariable under the 
reparametrization changes, this kind of transformations 
provokes the transition among macroscopic states as- 
sociated to different generalized canonical ensembles 
1281 which could non necessarily be mutually equivalent 
among them. 

The transformation rule EH allows us a better under- 
standing about the nature of the generalized thermostat 
associated to a generalized canonical ensemble by com- 
pared this last one with an ordinary thermostat dur- 
ing their interaction with a very large system. During 
a reparametrization change from a generalized canoni- 
cal ensemble with reparametrization (E) towards the 
Gibbs canonical ensemble the respective canonical pa- 
rameters n and (3 are related as follows: 



P(E; V ) 



90 (E ) 
dE 



-']■ 



(36) 



This equation tells us that a generalized thermostat with 
fixed canonical parameter n can be taken as an ordinary 



thermostat whose inverse temperature f3 depends on the 
instantaneous value of the total energy E of the system 
under analysis. Therefore, the total energy E and the 
inverse temperature (3 experience correlated fluctuations 
in the generalized canonical ensemble 1281 This fact is a 
remarkable feature of the generalized canonical ensemble 
in regard with the Gibbs canonical ensemble where p is 
fixed and E fluctuates around its mean value (E). We 
shall shown in the subsection IIII . El that the characteriza- 
tion of the fluctuations of the energy E and the inverse 
temperature [3 demands a suitable extension of the fluc- 
tuation theory of the conventional Thermodynamics to 
the present framework. 



D. Generalized Metropolis Monte Carlo method 

The basic problem in equilibrium statistical mechanics 
is to compute the phase space average, in which Monte 
Carlo method plays a very important role ifhl) . Among 
all admissible statistical ensembles used with the above 
purpose, the microcanonical ensemble provides the most 
complete characterization of a given system in thermody- 
namic equilibrium. However, microcanonical calculations 
could be difficult to carry out directly in a given applica- 
tion. A very simple way to overcome this difficulty is to 
consider the relationship among the microcanonical de- 
scription of a large enough system with other statistical 
ensembles (like the Gibbs canonical ensemble ^1 or the 
multicanonical ensemble |l5t fl6j) '). 

Suitable estimates of the microcanonical averages 
could be easily obtained from the Metropolis importance 
sample algorithm (MMC) if!?!) by using the equivalence 
between the microcanonical ensemble and the Gibbs 
canonical ensemble. This Monte Carlo method has some 
important features. It is extremely general and each 
moves involves O (1) operations. However, its dynam- 
ics suffers from critical slowing down, that is, if N is the 
system size, the correlation time r diverges as a criti- 
cal temperature is approached: the divergence follows a 
power law behavior r cx N z in second-order phase tran- 
sitions, while r diverges exponentially with iV in first- 
order phase transitions as consequence of the ensemble 
inequivalence |l8|) . 

The critical slowing down observed in the neighbor- 
hood of the critical point of the second-order phase tran- 
sitions is intimately related with the existence of the long- 
range order (divergence of correlation length) which char- 
acterizes the microscopic picture of this kind of phase 
transitions. This difficulty can be overcome in many 
physical systems by using appropriated clusters algo- 
rithms |l9|) . 

As already explained in the subsection IIII. Al the su- 
percritical slowing down (exponential divergence of the 
correlation times) of the first-order phase transitions is 
originated from the ensemble inequivalence of the canon- 
ical description. Such difficulty is usually overcome by 
using the multicanonical ensemble l)l5t Il6|) . a method- 
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ology which reduces the exponential divergence of the 
correlation times with respect to system size to a power 
in the first-order phase transitions. The multicanonical 
ensemble flattens out the energy distribution, which al- 
lows the computation of the density of states (E, N) 
for all values of E in only one run. 

An alternative way to overcome the supercritical slow- 
ing down of the MMC algorithm in the neighborhood of 
the first-order phase transitions is to ovoid the underlying 
ensemble inequivalence by using the generalized canoni- 
cal ensemble |2H1 The method, the generalized canonical 
Metropolis Monte Carlo (GCMMC) is explained in de- 
tails in the ref. so that, we shall limit in this subsec- 
tion to expose its more important features. 

Loosely speaking, the GCMMC algorithm is essen- 
tially the same MMC algorithm with the generalized 
canonical weight loqc (E) — exp [—77© (E)] instead of 
the usual Gibbs canonical weight cu c (E) = exp[—(3E]. 
The probability p for the acceptance of a Metropolis 
move is given now by p — min {1, exp [— 77 A0]} where 
A9 = Q (E + Ae) — 6 (E). Since Ae « E when the in- 
terest system is large enough, we can use the approxima- 
tion t/A© ~ {ridO (E) /dE} Ae = (3 (E; 77) Ae, which al- 
lows us to rewrite the acceptance probability p as follows: 

p~min{l,exp[-/3(E;r7) Ae]}. (37) 

The reader can recognize the natural appearance of the 
fluctuating inverse temperature (3 (E; 77) characterizing 
the generalized thermostat with parameter 77 and en- 
ergy reparametrization [E) introduced in the equation 
!5ol This feature of the GCMMC algorithm allows the 
Metropolis dynamics to explore regions which are inac- 
cessible for any other Monte Carlo method based on the 
Gibbs canonical ensemble, that is, regions which are char- 
acterized by exhibiting a negative heat capacity (see in 
subsection IIII.Af) . 

It is very easy to show why and when the GCMMC 
algorithm works. Considering the following expression of 
the generalized partition function: 

Z(V,N) = — I exp [-776 (E) + S (E, N)] dE, (38) 
Se J 

the application of the steepest descent method up to the 
Gaussian approximation allows us to express the con- 
ditions of the local ensemble equivalence between the 
macrostate of the generalized canonical ensemble 1281 with 
canonical parameter 77 and a microcanonical macrostate 
with total energy E e as follows: 

( SE T' = -j^2 s 'i " (£») - 5 " ( £ " w > > °- ( 40 > 

The condition 1391 identifies what it is usually re- 
ferred as the microcanonical inverse temperature (3 (E) = 



dS(E,N) /dE with the value of the effective inverse 
temperature (3 (E; 77) of the generalized thermostat at the 
stationary point E e . Reader may notice that this expres- 
sion is just the transformation rule 1351 The positive def- 
inition of the energy square dispersion (SE 2 ^ expressed 
in the equation I4UI imposes a restriction to the energy 
reparametrization (E) . It is very easy to show by com- 
bining the equations and 0U] that the positive defi- 
nition of the energy square dispersion (SE 2 ) is ensured 
by the concavity of the microcanonical entropy in the 
reparametrization 0: 

(SF*\- [ 9Q(E) \ 2 d 2 S ^ d 2 S 

We can used in principle several generalized canonical 
ensembles in order to describe all admissible macrostates 
within the microcanonical description. According to 
the ref. we can use the Gibbs canonical ensem- 
ble to describe all those energetic regions where the lo- 
cal ensemble equivalence takes place, while consider a 
generalized canonical ensemble with reparametrization 
Q(E) = Ntp(E/N): 

^^=exp{-A( £2 - £ )}, (42) 

within an anomalous region (e±, £2), with A -1 ~ /3 C , being 
e = E/N the energy per particle, and (3 C , an estimation 
of the inverse critical temperature of the first-order phase 
transitions within the Gibbs canonical ensemble. 



E. Thermodynamic uncertainly relation 

As already discussed, the total energy of the interest 
system E and the inverse temperature of the generalized 
thermostat (3 (E, 77) exhibit correlated fluctuations within 
the generalized canonical ensemble. Since the relative 
energy fluctuations 5E/E ~ 1/jV (with SE = ^(SE 2 )) 
drops to zero in the thermodynamic limit N — > 00, their 
corresponding microcanonical values in this limit are just 
the average expectation values within the generalized 
canonical ensemble E e = (E) and (3 (E e ) — ((3 (E, 77)). 

The correlation function (S(38E) between the fluctua- 
tions of the effective inverse temperature (3 (E; 77) of the 
generalized thermostat and the fluctuations of the total 
energy of the interest system SE can be obtained as fol- 
lows: 

5(3 ~ 770" (E e ) SE => (S/36E) = r?0" (E e ) (SE 2 ) , (43) 

which can be rephrased conveniently by eliminating the 
factor 77©" (E) using the equation The resulting ex- 
pression: 

{5{3SE) = l + (8E 2 ) d2Si f- N \ (44) 

is a very remarkable equation which could be referred as 
a thermodynamic uncertainly relation. 
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The reader can notice that the result^Jdoes not makes 
explicit reference to the reparametrization O (E), so that, 
it presents a general applicability for a large system in 
thermodynamic equilibrium interchanging energy with its 
surrounding. This identity represents a restriction among 
the fluctuations of the energy of the interest system E 
and the inverse temperature (3 of a given generalized ther- 
mostat, providing in this way a suitable generalization 
of the classical fluctuation theory lfl3h . Since the total 
energy E is just an admissible reparametrization for the 
microcanonical description, the reparametrization invari- 
ance supports without any lost of generality the following 
extension of the identity El 

(SnSe) = l + (Se 2 )^. (45) 

While effective inverse temperature of the general- 
ized thermostat can be fixed (6(3 — 0) whenever the 
entropy be a concave function in terms of the to- 
tal energy E, the inverse temperature dispersion 6(3 
can not vanish when the entropy be a convex func- 
tion. All those anomalous behaviors appearing in the 
Gibbs canonical ensemble during the ensemble inequiv- 
alence are also related with the downfall of the well- 
known Gaussian estimation of the energy dispersion 
(6E 2 ) = - {d 2 S(E,N) /dE 2 }' 1 in such anomalous re- 
gions. There d 2 S (E,N) /dE 2 > 0, and consequently, 
(5(3SE) > 1. A better analysis allows us to obtain 

the inequalities (SE 2 ) < - {d 2 S{E,N) /dE 2 }' 1 when- 
ever entropy be a concave function, while {6(3 2 ) > 
d 2 S (E, N) /dE 2 when entropy is convex. 

According to the identity^] the total energy E and 
the inverse temperature (3 are complementary thermody- 
namic quantities, in complete analogy like the time r and 
the energy E or the coordinates q and the correspond- 
ing momentum p are complementary quantities in the 
Quantum Mechanics (St5E > H, SqSp > H). Although 
there is nothing strange in the existence of this kind of 
relationship between temperature and the energy 2 , the 
conventional Thermodynamics usually deals with equi- 
librium situations where temperature is fixed and the 
energy fluctuates (Boltzmann-Gibbs ensemble), or the 
energy be fixed and the temperature fluctuates (micro- 
canonical ensemble). The generalized canonical ensemble 
1281 provides an appropriate framework for dealing with 
equilibrium situations where both energy and tempera- 
ture are not fixed. 

The thermodynamic uncertainly relation l44l admits the 
following generalization when the open system is con- 



2 Within the Quantum Mechanics framework the Gibbs canonical 
weight w (/3) = exp^ — flH^J can be considered as a evolution 

operator S (t) = exp I —irE/h) with imaginary time r —> —i/3h. 



trolled by several thermodynamical parameters: 

(6(3 k 6I m ) = S? + (6I n 6I m ) (46) 

expression with provides a new relevance of the entropy 
Hessian within the present generalization of the classical 
fluctuation theory. The development of the present ideas 
deserves a further study. 

F. The role of the external conditions 

The reader can noticed by a simple inspection of the 
stationary conditions [M] that the local ensemble equiva- 
lence demands the local concavity of the microcanonical 
entropy, while the global ensemble equivalence demands 
the global concave character of this thermodynamic po- 
tential. As already illustrated in the subsections IIII.BI 
and IIII. CI such conditions depend crucially on the na- 
ture of the reparametrization (E) used in the general- 
ized canonical ensemble 1281 

Obviously, regions of ensemble inequivalence also 
depends on the reparametrization. Since different 
reparametrizations in the generalized canonical ensem- 
ble 1281 represent different experimental setups (the us- 
ing of different generalized thermostats), the existence of 
all those anomalous behaviors associated to the ensem- 
ble inequivalence depend crucially on the nature of the 
experimental setup. Therefore, we can not refer in the 
present framework to the existence of an ensemble equiv- 
alence without indicating the external conditions which 
have been imposed to the interest system. This situation 
is completely analogue to the question about the charac- 
ter of the motion of a given particle without specifying 
the reference frame. 

As already pointed in the subsection IIII. Al the phe- 
nomenon of phase coexistence associated to the ensemble 
inequivalence is interpreted as a first-order phase tran- 
sition in the conventional Thermodynamics. However, 
the ensemble inequivalence represents an anomaly within 
the canonical description which reflects the inability of 
this ensemble in considering all those macrostates ob- 
served in the microcanonical description of the interest 
system. Consequently, what we called a first- order phase 
transition in the conventional Thermodynamics is a phe- 
nomenon only relevant in the Gibbs canonical descrip- 
tion: such behavior does not represent anything anoma- 
lous within the microcanonical description. 

Obviously, the presence of an ensemble inequivalence 
in the Gibbs canonical description can be easily recog- 
nized within the microcanonical ensemble by the exis- 
tence of convex regions of the entropy when this de- 
scription is performed in terms of the total energy E. 
However, it is necessary to emphasize that there is noth- 
ing anomalous in this behavior of the entropy since the 
microcanonical description is reparametrization invari- 
ance, and consequently, any anomaly within this ensem- 
ble must be also reparametrization invariant. Therefore, 
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the convex or concave character of the entropy can not be 
used as an indicator of real anomaly within the micro- 
canonical description since this mathematical property 
can be arbitrarily modified with a simple change of the 
reparametrization E — > O (E) . 

It is convenient to recall one time again the generic (ex- 
perimental) definition of phase transition: " ... a phase 
transition brings about a sudden change of the macro- 
scopic properties of a system while smoothly varying of a 
control parameter" . To be more precise, we need to add 
that the behavior of the system is affected by the way 
its thermodynamic equilibrium is controlled by the ex- 
perimental setup, and consequently, the thermodynamic 
anomalies which could be appearing under certain ex- 
ternal conditions may not observed by using other ex- 
perimental setup. Furthermore, taking into account the 
example of the Gibbs canonical ensemble which is un- 
able to describe many thermodynamic states during the 
occurrence of a first order phase transition, it is very 
important to remark that other relevant thermodynamic 
anomalies could be hidden behind of any lost of informa- 
tion associated to the ensemble inequivalence. 



IV. A PARADIGMATIC EXAMPLE 

A. Hamiltonian of the Potts model 

The present section will be devoted to perform the 
thermostatistical description of a simple model exhibit- 
ing phase transitions in order to illustrate some aspects 
already discussed in the above section. Let us be gan the 
present section by introducing the Potts model (|2Lr) . a 
toy model whose Hamiltonian is defined as: 

^ = EI 1 -W, ( 47 ) 

on a hypercubic dimensional lattice, being Oi the spin 
state at the i-th lattice point which can take q (> 2) pos- 
sible integer values or components, <Tj = 1, 2, ...q, where 
the sum is over pairs of nearest neighbor lattice points. 
Hereafter we will assume periodic boundary conditions. 
It is very easy to notice that the Potts model with q = 2 
is equivalent to the Ising model of the ferromagnetism. 
Obviously, the cases with q > 2 are suitable generaliza- 
tions of this paradigmatic toy model which also admit a 
ferromagnetic interpretation (see in subsection II V . C jl . 

Potts models are very amenable for studying phase 
transitions with different order, i.e. q = 2 (Ising model) 
exhibits a continuous (second-order) phase transition, 
while q = 10 shows a discontinuous (first-order) phase 
transition. The present section will be devoted to study 
the main thermodynamic features of the q = 10 Potts 
model on a L x L square lattice by using the GCMMC 
algorithm described in the subsection IIII.DI and the 
ref.©. 



B. Thermodynamical potentials for q = 10 

The q — 10 states Potts model exhibits a first-order 
phase transition associated to the ensemble inequiva- 
lence, which provokes the existence of a supercritical 
slowing down during the ordinary Metropolis dynamics 
based on the consideration of the Gibbs canonical ensem- 
ble. Clusters algorithms, like Swendsen-Wang or Wolf 
algorithms (19), do not help in this case due to they are 
still based on the consideration of the canonical ensem- 
ble 1)221) . and consequently, the supercritical slowin g do wn 
associated to the ensemble inequivalence persists (fill) . 

A preliminary calculation by using the MMC algorithm 
allows us to set a anomalous region within the energetic 
windows (61,62) with 6\ = 0.2 and 62 = 1-2 in which 
takes place the critical slowing down associated with the 
existence of a first-order phase transition in this model 
for L = 25. The inverse critical temperature was esti- 
mated as [3 C ~ 1.4, allowing us to set A = 0.8. The 
main microcanonical observables obtained by using the 
GCMMC algorithm with n — 10 5 Metropolis iterations 
for each points is shown in the FIG^ th e entropy per 
particle s (e), the caloric (3 (e) = ds (e) /cte and curvature 
k (e) = d 2 s(6)/de 2 curves. A comparative study be- 
tween the present method and the Swendsen-Wang clus- 
ters algorithm iflU) is shown in the FIG [3 

Notice that the heat capacity c (e) = — 1 (e) / n (e) 
becomes negative when e £ (e Q ,£fc) with e a ~ 0.51 and 
6b — 0.93, which is a feature of a first-order phase tran- 
sition in a small system becoming extensive in the ther- 
modynamic limit. Such anomalous behavior is directly 
related with the backbending in the caloric curve (3 versus 
6 and the existence of a convex intruder in the relative 
microcanonical entropy per particle As versus 6 shown 
in the inserted graph. This last dependence was obtained 
from a simple numerical integration of the caloric (3 (e) 
and curvature k (e) dependences by using the scheme: 

Ss (e) = s (e + fe) — s (e) , 

~ (3 (6)66 + ^k(6)S6 2 , (48) 

based on the second-order approximation of the power 
expansion of the entropy. We clarify to the reader that 
the true dependence plotted in this figure is given by 
As* (e) = s* (e) - s* (e*), where s* (e) = s (e) - ae, with 
£* = 0.015 and a = 1.41, in order to make more evidence 
the existence of the convex intruder. It can be shown that 
the convex intruder disappears progressively with the in- 
creasing of the system size until becoming the Maxwell 
line represented in the inserted graph of the FIG^ which 
provides us information about the transition temperature 
and the latent heat of the first-order phase transition un- 
dergone by this model system. 

Since the heat capacity is always positive within the 
canonical ensemble, such anomalous regions are inacces- 
sible in this description. This fact evidences the existence 
of a significant lost of information about the thermo- 
dynamical features of the system during the occurrence 
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FIG. 1 Microcanonical caloric curve (5 (e) = ds (e) /de and 
curvature k (e) = d 2 s (e) /de 2 of the 25 x 25 Potts model with 
q = 10 states with periodic boundary conditions obtained 
by using the GCMMC algorithm. Notice the energetic re- 
gion with a negative heat capacity. The inserted graph shows 
the behavior of the relative microcanonical entropy As (e) in 
which is remarkable the existence of a convex intruder. Errors 
are smaller than the symbols linear dimension. 



of a first-order phase transitions when the canonical en- 
semble is used instead of the microcanonical one. This 
difficulty is successfully overcome by the GCMMC algo- 
rithm, which is able to predict the microcanonical av- 
erage of the microscopic observables in these anomalous 
regions where any others Monte Carlo methods based on 
the consideration of the Gibbs canonical ensemble such 
as the original MMC, the Swendsen-Wang and the Wolff 
single cluster algorithms certainly do not work. 

This fact is clearly illustrated in the FIG0 which 
evidences that the Swendsen-Wang algorithm is unable 
to describe the thermodynamic states with a negative 
heat capacity: its results within the anomalous region 
differ significantly from the ones obtained by using the 
GCMMC algorithm. The Swendsen-Wang dynamics ex- 
hibits here an erratic behavior originated from the com- 
petition of the two metastable states present in the 
neighborhood of the critical point (the energy distribu- 
tion function in the canonical ensemble is bimodal when 
P G (/3i,/3 2 ), where Pi = 1.405 and (3 2 = 1.445, being 
this feature the origin of the supercritical slowing down) . 

As already shown by Gross in ref . Qj , a convex intruder 
can be associated with the existence of a non-vanishing 
interphase surface tension during the first-order phase 
transitions in systems with short-range interactions out- 




FIG. 2 Comparative study between the GCMMC and the 
Swendsen-Wang (SW) cluster algorithm. The SW algorithm 
is unable to describe the microcanonical states with a negative 
heat capacity. Notice the bimodal character of the energy 
distribution function at f3 = 1.42. 



side the thermodynamic limit. Nevertheless, we will show 
in the next subsection that this is not the only one infor- 
mation which could be hidden behind the presence of a 
negative heat capacity. 

We show in FIG^l the canonical description of this 
model: the caloric curve and the microcanonical Planck 
thermodynamic potential per particle p m (e) versus p (e) 
obtained from the Legendre transformation p m (e) = 
(3 (e) e — s (e). The branches OA and AO' represent the 
canonically stable thermodynamic states corresponding 
to the minimal values of the microcanonical Planck po- 
tential at a given /3 (the true Planck thermodynamic po- 
tential per particle of the canonical description when the 
system size N is large enough can be approximated by: 



Pc (p)=-lnZ(p,N) /N 



min \p ri 



(e)}, (49) 



being Z (P,N) the canonical partition function). Thus, 
the points A represent the critical point of the first-order 
phase transitions where Pa — 1-421, which is recognized 
from the condition: 

P A = (3 (etf) = P and p (e { A ] ) = p , 

and appreciated in the inserted graph. Obviously, the 
first derivative of the canonical Planck thermodynamic 
potential 1491 exhibits a discontinuity at P — Pa- Since 



-(-) 



e A ' ~ 1.099 and e A 



(+) 



0.319, the estimated latent heat 



qih associated to this phase-transition is given by qih 



e ( A ) -e ( ~t ] - U.7N. 



The energy distribution function within the canoni- 
cal ensemble is bimodal inside the interval {Pc , Pb) with 
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FIG. 3 Canonical description of the q — 10 states Potts model 
on 25 x 25 square lattice in which are represented the caloric 
curve e versus j3 and the Planck thermodynamic potential 
per particles P versus /3 dependences. The inserted graph is 
a magnification of the region enclosed in the circle indicating 
the exact point of the first-order phase transition. 



(3b — 1-405 and (3c — 1-443, which is directly related 
with the existence of metastable states, the branches AC 
(supercooled disordered states) and AB (superheated or- 
dered states), which are the origin of the supercritical 
slowing down behavior observed during the Monte Carlo 
simulations by using the ordinary Metropolis importance 
sample algorithm MMC (the exponential divergence of 
the correlation time with the system size increasing). 

The branch BC is canonically unstable since its points 
represent thermodynamical states with a negative heat 
capacity associated to the convex intruder of the mi- 
crocanonical entropy, where es — 0.933 and eq — 
0.512. The energy region (ec,£b) is practically invisi- 
ble within the canonical description when the system is 
large enough, and consequently, the ordinary Metropo- 
lis algorithm based on the Gibbs ensemble will never 
access there as a consequence of the ensemble inequiv- 
alence. Fortunately, the Metropolis algorithm based on 
the reparametrization invariance (the using of general- 
ized canonical ensembles) considered in the present study 
overcomes successfully all those the difficulties undergone 
by the original method 



C. Magnetic properties for q = 10 

The Potts model can be easily rephrased as ferromag- 
netic system with the introduction of the bidimensional 



vector variable Sj = [cos(K<7i) , sin(/ccrj)] with k = 2ir/q, 
where the total magnetization is simply given by M = 
Y]. Sj. Their magnetic thermostatistical properties will 
be studied in the present subsection by means of the 
average magnetization (M) and the square dispersion 

G = ^(M — (M)) 2 ^} which characterizes the fluctuations 

of this microscopic observable. 

It is well-known that the square dispersion of the to- 
tal magnetization G provides us a measure g about the 
spatial correlations between the spin particles 



G 
N 



N ^ 



9ij: 



(51) 



where g L j — (sj • Sj) — (s,) • (sj) is the two-point corre- 
lation function. This quantity is also related with the 
magnetic susceptivity x = dM/dB of the system under 
the presence of an external magnetic field B (coupled to 
the total magnetization modifying the original Hamilto- 
nian H S3 as follows H B = H - B • M) throughout the 
called static susceptibility sum rule (0): 



dP_ 
~dB 



(3M, 



dM 
~dB 



PG C , 



(52) 



within the canonical description, while the corresponding 
identities in the microcanonical ensemble are given by: 



dS 01# dM a „ 
- = /3M ! — = ^ 



M- (53) 



where P and Sb are the Planck potential of the canon- 
ical ensemble and the Boltzmann entropy of the micro- 
canonical ensemble respectively, B = |B| and M is the 
projection of the total magnetization along the external 
magnetic field vector (see demonstration in appendix[nj|. 

It is very important to remark that although the exis- 
tence of the ensemble equivalence in the thermodynamic 
limit allows us to identify asymptotically the expecta- 
tion values of the microscopic observables (like the total 
magnetization M), those quantities characterizing their 
fluctuations (G) and the response functions (%) depend 
essentially on the nature of the statistical ensemble used 
in the description. For example, while the canonical sus- 
ceptibility is always positive, Xc = PG C > 0, its corre- 
sponding microcanonical quantity Xm could be negative 
as a consequence of the presence of the term MdM/dE 
in the second relation of the equation 1531 since magneti- 
zation decreases in the ferromagnetic systems with the 
energy increasing. 

Let us return to the study of the magnetic proper- 
ties of the Potts model. It is very easy to see that the 
Hamiltonian 1471 is invariant under the discrete group of 
transformations LT g composed by the q\ permutations 
among q spin states. The subgroup of cyclic permu- 
tations Z q — {Tfe; k = 1,2, ...q} is also a subgroup of 
the group of bidimensional rotations U (2), whose lo- 
th transformation acting on a given vector variable 
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FIG. 4 Magnetization m versus j3 dependence o the q = 10 
states Potts model with L — 25, which shows the existence of 
a discontinuous ferromagnetic-paramagnetic phase transition 
at p = Pa = 1.421. 



s (cr) = [cos (kit) , sin (ko~)] induces a rotation in an angle 
Aip — 2irk/q: 



T k s (a) = s(a + k) 



(54) 



Obviously, the existence of this last symmetry leads triv- 
ially to the vanishing of the exact statistical average of 
the magnetization. 

However, a net magnetization can appears at low ener- 
gies during the paramagnetic-ferromagnetic phase tran- 
sition (from disordered states towards the ordered ones) 
as a consequence of the spontaneous symmetry break- 
ing associated to the nontrivial occurrence of an ergodic- 
ity breaking in the underlying dynamical behavior of this 
model system. This dynamical phenomenon manifests it- 
self when the time averages and the ensemble averages of 
certain macroscopic observables can not be identified due 
to the microscopic dynamics is effectively trapped in dif- 
ferent subsets of the configurational or phase space dur- 
ing the imposition of the thermodynamic limit N — > oo 



We show in the FIG0]the magnetization density m = 
|(M)| /N within the canonical ensemble, where the sta- 
ble (OA and AO') and metastable branches (AC and 
AB) are clearly visible. The existence of a ferromagnetic- 
paramagnetic phase transition at the critical point (3a 
shows clearly that the Z q symmetry has been sponta- 
neously broken. However, it is remarkable how a net 
magnetization appears abruptly, that is, with a discon- 
tinuous character, a behavior which is very similar to 
the one observed during the solid-liquid first-order phase 
transition where the traslational symmetry is also spon- 
taneously broken. 

As already commented in the above subsection, there 
are many thermodynamic information which could be 
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FIG. 5 Magnetic properties of the model within the micro- 
canonical ensemble in which is shown the existence of two 
phase transitions at eff — 0.7 (ferro-ferro) and £/ P — 0.8 
(ferro-para). 



hidden behind a negative heat capacity during the canon- 
ical description. Particularly, the reader may agree with 
us that it is reasonable to expect within the microcanon- 
ical description the existence of a critical energy where 
ferromagnetic-paramagnetic phase transition takes place 
without the discontinuous character of the magnetization 
curve observed within the canonical description. Surpris- 
ingly, reader can notice that the qualitative behavior of 
the magnetic properties of this model system shown in 
the FIG|2l (magnetization density m and the dispersion 
g = G/N) are much more interesting than our prelimi- 
nary idea. 

This magnetization density m versus e dependence ev- 
idence clearly what could be considered as two phase 
transitions within the microcanonical description of the 
q = 10 states Potts model for L = 25: a continu- 
ous (para-ferro) phase transition at the critical point 
£f P — 0.8, and a discontinuous (ferro-ferro) phase tran- 
sition at Eff ~ 0.7. Most of thermodynamic points were 
obtained from a data of n = 10 5 Metropolis iterations, 
with the exception of all those points belonging to the 
energetic interval (0.7,0.93) where n = 5 x 10 6 iterations 
were needed in order to reduce the significant disper- 
sion of the expectation values close to the critical points. 
The large dispersions observed throughout the disper- 
sion g, the large relaxation times during the Metropolis 
dynamics, and the qualitative behavior of the magneti- 
zation density strongly suggest us the presence of several 
metastable states with different magnetization states at 
a given energy within this last region, which is demon- 
strated in the FIG0 

The microcanonical continuous phase transition be- 
tween the paramagnetic-ferromagnetic type I (with low 
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FIG. 6 Histogram of the projection of the microscopic magne- 
tization along the direction of the spontaneous magnetization 
of the ferromanet type I, which demostrates the existence of 
metastable states when e = 0.844. 



magnetization density) phases takes place with the spon- 
taneous breaking of the Z q symmetry associated to the 
occurrence of an ergodicity breaking in the microscopic 
dynamics. Apparently, the large fluctuations and long- 
range correlations ordinarily associated to this kind of 
phase transition are overlapped with the large fluctua- 
tions associated to the existence of metastable states. 

On the other hand, the Z q symmetry has been already 
spontaneously broken during the occurrence of the mi- 
crocanonical discontinuous phase transition between the 
ferromagnet type I - ferromagnet type II (with high 
magnetization density). However, this phase transition is 
also associated to the occurrence of an ergodicity break- 
ing in the underlying microscopic picture: the micro- 
scopic dynamics of this model system can be effectively 
trapped in the thermodynamic limit L — ► oo in any of 
the metastable states with different magnetization den- 
sity present in the neighborhood of this critical point. 

The FIGs[7| and [S] show the magnetic properties of the 
model system under the presence of an external magnetic 
field with B = 0.01 and B = 0.1. The reader can no- 
tice in the FIG[7|that the microcanonical discontinuous 
phase transition is also present for low intensities of the 
external magnetic field, but the discontinuity observed in 
the magnetization m (e) and the dispersion g (e) depen- 
dences is not so abrupt in this case. Apparently, such 
behavior is reduced progressively with the B increasing 
until disappear when the magnetic field intensity is large 
enough, as already illustrated in the FIG 00 This last 
figure shows that the dispersion curve has a peak, which 
is a remanence of the paramagnetic-ferromagnetic phase 
transition when B 0. Reader may notice that the 
presence of discontinuous phase transition affects signif- 
icantly the qualitative behavior of the square dispersion 
g of the magnetization. 
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FIG. 7 Magnetic properties of the model under the pres- 
ence of an external magnetic field with B = 0.01. Notice 
that in this case the dispersion is not large than in the first 
case with B = 0, but it is clearly visible the existence of the 
ferromagnetic-paramagnetic phase transition. 
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FIG. 8 The same dependences of the FIGs |0 and 0, but 
now under the influence of an external magnetig field with 
B = 0.1. 



D. What can be learned from the microcanonical 
description of this model system? 

The backbending behavior in the microcanonical 
caloric curve shown in the FIG^ (becoming a plateau in 
the thermodynamic limit) is usually interpreted as the 
signature of a first-order phase transition (with spon- 
taneous symmetry breaking). On the other hand, the 



14 



magnetization curve shown in the FIG [5] exhibits two 
anomalies which can be also considered as phase tran- 
sition of the thermostatistical description of this model 
system: a continuous phase transition (with spontaneous 
symmetry breaking) and a discontinuous phase transi- 
tion (where the underlying symmetry has been already 
broken). Consequently, what kind of phase transition 
exhibits the thermostatistical description of the q = 10 
states Potts model? The answer of the above question in 
our viewpoint depends on the external conditions imposed 
to the system. 

Obviously, this model system put in contact with a 
Gibbs thermostat undergoes a discontinuous first order 
phase transition while a smooth change of the inverse 
temperature (3 in the neighborhood of the critical point 
Pa = 1.421. There, the phase coexistence phenomenon 
characteristic of this kind of phase transitions can be ob- 
served and the region (ec,£s) is inaccessible within the 
Gibbs canonical description in the thermodynamic limit. 

A different picture is revealed when the model system 
is isolated (microcanonical description). All the anoma- 
lous region (ec, £b) is now accessible. Many microcanon- 
ical states there show no anomalous behavior with the ex- 
ception of those macrostates within the region (g^*, Eg) 
which are affected by the incidence of metastable states. 
The imposition of the thermodynamic limit leads to the 
suppression of the metastable states, but the existence 
of anomalous behaviors persists in the neighborhood of 
the para-ferro continuous phase transition at e f p and the 
ferro-ferro discontinuous phase transition at £//. Inter- 
estingly, such anomalies can not be apparently associated 
with any anomaly of the caloric or curvature curves in 
the FIG^ but by using the microcanonical magnetiza- 
tion curve shown in FIG [5] 

Taking into account the microcanonical thermody- 
namic identities in the thermodynamic limit N — > oo 
shown in the equation 1531 between the entropy S, the 
magnetization M, the external magnetic field B and the 
square dispersion G: the discontinuous phase transition 
at Sff corresponds to a discontinuity of the first deriva- 
tive of the entropy per particle: 



r r ds{e,B) ds{e,B) 
nm nm — i= urn lim — 

^s+B^O dB ^ e-^jB^O dB 



(55) 



while the continuous phase transition at e f v corresponds 
to a discontinuity of the second derivative of the entropy: 



hm lmi — - ^ lim lim — - . (56) 

,^ £ + b^o dedB e ^ e - b^o dedB 

' s fp e e ip 



Consequently, the present microcanonical anomalies can 
be recognized by the lost of analyticity of the entropy per 
particle in the thermodynamic limit, and the same ones 
are related with the occurrence of an ergodicity breaking 
in the microscopic picture of this model system. 



V. TOWARDS A NEW CLASSIFICATION 
SCHEME 

The analysis developed in the above two sections al- 
lows us to classify the thermodynamic anomalies type A 
and type B introduced in the beginning of the section 
IIIII as follows. The anomaly type A are just anomalies 
observed in open systems (canonical description) which 
are dependent on the external experimental conditions 
due to they originated from the ensemble inequivalence. 
The anomaly type B are all those anomalies which are 
always observed in an isolated system (microcanonical 
description) and can potentially appear under any exter- 
nal experimental conditions. 



A. The anomalies type A 

Anomalies type A correspond what we usually called 
as first-order phase transitions in conventional Thermo- 
dynamics, and they are only relevant for an open system 
since their existence depends crucially on the nature of 
the external conditions imposed to the interest system. 
We recognized them by the existence of a latent heat for 
the phase transition associated to the multimodal char- 
acter (a signature of the phase coexistence) of the en- 
ergy distribution function within the (generalized) canon- 
ical description. In terms of the thermodynamic poten- 
tials, an anomaly type A manifests as a discontinuity in 
some of the first derivatives the generalized Planck po- 
tential, which are directly related to the existence of con- 
vex regions of the microcanonical entropy S in a given 
reparametrization O dependent on the external experi- 
mental setup (the using of a generalized thermostat). 

An anomaly type A represents in this way the inabil- 
ity of the (generalized) canonical description in describ- 
ing all those macrostates which can be accessed within 
the microcanonical description. As already shown, this 
kind of anomaly can be avoided by using another experi- 
mental setup which ensures the global equivalence of the 
generalized canonical description with the microcanoni- 
cal ensemble. The Monte Carlo method exposed in the 
subsection IIII.DI is precisely based on this idea. 

The Swensen-Wang algorithm is a Monte Carlo clus- 
ter algorithm based on the consideration of the Gibbs 
canonical ensemble, that is, it simulates the thcrmody- 
namical equilibrium of the Potts model system put in 
contact with an ordinary thermostat (heat bath). Since 
the Gibbs canonical description of the q = 10 states Potts 
model is not globally equivalent to its microcanonical de- 
scription, a sudden change in the energy per particle e of 
this model system is observed in the neighborhood of the 
critical inverse temperature /3 C ~ 1.421. On the other 
hand, the GCMMC algorithm describes the thermody- 
namic equilibrium of a given system in contact with a 
generalized thermostat. As already shown in FIG El this 
change of the external conditions allows the system to 
access to all those inaccessible regions within the Gibbs 
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FIG. 9 Schematic representation of the onset of the magneti- 
zation curve within the microcanonical description in the limit 
of zero applied magnetic field of a ferro-para phase transition. 



canonical ensemble, and therefore, neither there exist lost 
of information nor any sudden change of the energy per 
particle is observed now. Thus, the thermodynamical 
study of the q — 10 states Potts model shows that the 
first-order phase transitions are avoidable thermodynam- 
ical anomalies. 



B. The anomalies type B 

The lost of analyticity of the entropy per particle in the 
thermodynamic limit, the anomaly type B, is obviously 
the only one mathematical anomaly of the microcanon- 
ical entropy which is reparametrization invariant. This 
character explains why this kind of anomaly could po- 
tentially appear under any external experimental setup 
(which determines the specific reparametrization used 
in the generalized canonical ensemble), and consequently, 
it is the only thermodynamical anomaly microcanonically 
relevant. We say "potentially" because of there exists 
the possibility that an anomaly type B could be hid- 
den by the lost of information associated to the ensemble 
inequivalence, i.e.: the continuous (ferro-para) and the 
discontinuous (ferro-ferro) phase transitions described in 
the FIG 03 are hidden within the Gibbs canonical descrip- 
tion. 

Since the microcanonical ensemble is just a dynami- 
cal ensemble, the anomaly type B should be the macro- 
scopic manifestation of a sudden change in the dynami- 
cal behavior of the isolated Hamiltonian system. As al- 
ready pointed in our study of the q = 10 states Potts 
model, the lost of analyticity of the entropy per particle 
in the thermodynamic limit seems to be related to the 
ergodicity breaking phenomenon. The well-known spon- 
taneous symmetry breaking associated to the ferro-para 
second order phase transitions schematically represented 
in the FIG 03 is a generic example of an ergodicity break- 
ing which is always connected to a lost of analyticity of 
the entropy per particle in the thermodynamic limit. The 
demonstration starts from the consideration of the first 



ds (e, B) 
dB 



(3{s,B)m(e,B), 



(57) 



which allows us to obtain the magnetization curve 
m (e, B) from the entropy per particle in the thermo- 
dynamic limit s (e, B), being (3 (e, B) = ds (e, B) jde the 
microcanonical caloric curve. The partial derivative d/de 
of the above equation is given by: 



d 2 s (£, B) _ d(3 (e, B) 



dedB 



ds 



m [e, B)+P (e, B) 



dm [e, B) 



de ■ 

(58) 

While the first derivative of the caloric curve 
d/3 [s, B) jde always exists in short-range interacting sys- 
tems, the dependence dm [s, B) jde exhibits a disconti- 
nuity at the critical point s c in the limit of zero applied 
magnetic field, and consequently: 



Me)=hm^M 
v ; b-+o dedB 



(59) 



is discontinuous function. 

The anomaly type B described above can be referred 
as a microcanonical continuous phase transition since the 
first derivatives of the entropy per particle in the thermo- 
dynamic limit are continuous [Anomaly type B.I). Con- 
trary, we can referred a anomaly type B as a microcanon- 
ical discontinuous phase transition when some of the first 
derivatives of the entropy per particle in the thermody- 
namic limit are discontinuous at the point of lost of ana- 
lyticity [Anomaly type B.II). Particularly, the ferro-ferro 
phase transition observed in the microcanonical descrip- 
tion of the q = 10 states Potts model is a clear example 
of a microcanonical discontinuous phase transition: Since 
the magnetization curve is discontinuous, and therefore, 
the first derivative 1571 is also discontinuous. 

The reader can notice that most of the well-known 
continuous (second-order) phase transitions in the con- 
ventional Thermodynamics corresponds to microcanoni- 
cal continuous phase transitions, since the applicability 
of the Legendre transformation P = [3E — S during the 
ensemble equivalence in the thermodynamic limit allows 
that every anomaly type B.I leads to lost of analyticity 
of the Planck potential (or the Helmholtz free energy). 
However, we will show below that all anomaly that could 
be classified as a continuous phase transition within the 
conventional Thermodynamics does not correspond nec- 
essarily to a microcanonical continuous phase transition. 

A feature of the physical systems exhibiting a contin- 
uous phase transition associated to the occurrence of an 
spontaneous symmetry breaking is the existence of diver- 
gent power laws behavior in the heat capacity (and other 
response functions) in the neighborhood of the critical 
point C : 



c(f3) 




(60) 
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being t = j3 c — j3, where the called critical exponent a 
and the amplitude ratio /A^ are the same within 
a universality class (0). The characteristic "A-form" of 
the continuous phase transitions is directly related to 
the presence of divergent power laws with universal ratio 
A^ + ' /j4>) 7^ 1 and nontrivial critical exponent a: i.e.: 
a = (0.110 — 0.116) for the Ising universality class. The 
divergence of the heat capacity in the continuous phase 
transitions can be explained by an eventual vanishing of 
the second derivative of the entropy d 2 s/de 2 at the crit- 
ical point. However, the nature of the divergent power 
law depends crucially on the analyticity of the entropy 
per particle in the thermodynamic limit. 

Let be a hypothetical system whose entropy per par- 
ticle in the thermodynamic limit s (s) is analytical in a 
given region but its second derivative d 2 s (s) /ds 2 van- 
ishes eventually at a certain point e c of this region. This 
is just anomaly type A where the anomalous region of the 
microcanonical entropy is composed by only one point. 
Since the caloric curve /3(e) = ds (e) /de is bijective in 
this case, ensemble equivalence is ensured. However, the 
average square dispersion of the system energy and the 
heat capacity within the Gibbs canonical ensemble go to 
infinite at e c , so that, the second derivative of the Planck 
thermodynamic potential diverges at the corresponding 
critical point f3 c = (3 (e c ). In spite of the present anoma- 
lous behavior can be classified as a continuous phase tran- 
sition in the conventional Thermodynamics viewpoint, it 
does not correspond to any microcanonical phase transi- 
tion. Particularly, the missing of the lost of analyticity 
of the entropy allows us to think that the above anomaly 
can not be associated to the occurrence of an ergodic- 
ity breaking. For example, the analytical character of 
the entropy per particle in the thermodynamic limit of 
the above hypothetical system allows us to approximate 
the caloric curve (3 (e) in the neighborhood of the critical 
point e c by the Taylor power series: 

P = f3 c -\(e-e c ) 2n+1 + ..., (61) 

where A > and n a positive integer, expression lead- 
ing to a power law divergent behavior of the heat ca- 
pacity E01 with critical exponent a — 2n/(2n+l) and 
A(~) = A (+ ) = /?o/A 1_Q (2n + l). Consequently, this re- 
sult corresponds to a critical phenomenon with universal 
ratio A^ /A^ = 1 and critical exponent a related to a 
odd number In + 1 = 1/ (1 — a). 

The above result allows to claim that the existence of 
divergent power laws with nontrivial critical exponents a 
and universal ratio A^ /A^ =/= 1 in most of real phys- 
ical systems exhibiting a continuous phase transition as- 
sociated to the occurrence of a spontaneous symmetry 
breaking is a clear indicator of the lost of analyticity of 
the entropy per particle in the thermodynamic limit in 
such cases. An anomaly like the one exhibited by our hy- 
pothetical system can be seen as an asymptotic case of 
a discontinuous phase transition, which shall be referred 
in this work as a Ist-kind continuous phase transition. 

The discontinuous character of the first derivatives of 



the entropy per particle in the thermodynamic limit dur- 
ing the occurrence of the microcanonical discontinuous 
phase transitions can lead to discontinuity of some of the 
first derivatives of the Planck potential (or other ther- 
modynamic potential characterizing an open system) or 
even provoke the discontinuity of the thermodynamical 
potential itself. This last mathematical anomaly is un- 
usual for the systems dealt within the conventional Ther- 
modynamics, but it can be observed in the astrophysical 
context and others long-range interacting systems. The 
interested reader can find in the ref. \2"^ \ an example of 
a discontinuity of the caloric curve [3 (e) observed in the 
thermodynamical description of astrophysical model, in- 
dicating the existence of metastable states at the same 
total energy with different temperature (mathematical 
anomaly leading to a discontinuity of the Planck poten- 
tial P = (3E — S) . While the thermodynamical behavior 
associated to the discontinuity of the Planck thermody- 
namic potential could be referred as a zero-order phase 
transition within the well-known Ehrcnfcst classification, 
the discontinuity of some of its first derivatives is just a 
discontinuous phase transition. 

Although we can not provide in this work a rigorous 
demonstration about the relationship between the lost 
of analyticity of the entropy per particle in the thermo- 
dynamic limit with the occurrence of ergodicity break- 
ing in the microscopic dynamics, we have considered in 
this work some examples suggesting that such connec- 
tion exists. Loosely speaking, the ergodicity breaking 
takes place as a consequence of the dynamical competi- 
tion among metastable states, where the predominance of 
any of them depends crucially on the initial conditions of 
the microscopic dynamics and the boundary conditions 
(0; 0). Particularly, the thermodynamical study of the 
q = 10 states Potts model presented in the section llVl sug- 
gests that the microcanonical continuous phase transition 
observed in this model system is provoked by the compe- 
tition among metastable states which essentially identical 
(states with the same magnetization density) since they 
are related by a symmetry transformation (Z q ). On the 
other hand, the microcanonical discontinuous phase tran- 
sition observed in this model system is also provoked by 
the competition among metastable states, which are es- 
sentially different (states with different spontaneous mag- 
netization density) since them can not be related by a 
symmetry transformation. 

The reader may notice that the discontinuity of the 
first derivative of the entropy described in the ref. l|23|) 
(a microcanonical discontinuous phase transition) can be 
associated to the occurrence of ergodicity breaking orig- 
inated from the competition among metastable states 
with different temperature. 



C. Summary 

As already shown in this work, the existence and the 
features of anomalies observed in the thermodynamical 
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Thermodynamical anomaly 


Classification 


Type A (presence of non 
concave regions of the entropy S) 
(non reparametrization invariant) 


microcanonically 
irrelevant 


Type B.I (lost of analyticity of 
the entropy S with continuous 
first derivatives) 
(reparametrization invariant) 


microcanonical 
continuous PT 


Type B.II (lost of analyticity of 
the entropy S with some 
discontinuous first derivatives) 
(reparametrization invariant) 


microcanonical 
discontinuous PT 



TABLE I Phase transitions (PT) within the microcanonical 
description (isolated system). 



description of given Hamiltonian system depends cru- 
cially on the external conditions which have been im- 
posed. It means that a phase transition is not only an in- 
trinsic thermodynamical anomaly of a given system, but 
also a specific response to the nature of the external con- 
trol of its thermodynamic equilibrium. Generally speak- 
ing, the mathematical description of the phase transitions 
for systems in thermodynamic limit starts from the con- 
sideration of the lost of analyticity of the thermodynamic 
potential which is relevant in a given application. 

The most simple characterization of the thermodynam- 
ical anomalies of a given Hamiltonian system is obtained 
within the microcanonical description, which is relevant 
when the interest system is isolated. Phase transitions 
here can be recognized by the lost of analyticity of the 
entropy S in the thermodynamic limit, a mathematical 
anomaly (type B.I or B.II) which is reparametrization 
invariant and should be originated from the occurrence 
of the ergodicity breaking phenomenon. This results are 
summarized in the Table |U 

The presence of the experimental setup controlling 
the thermodynamic equilibrium of the interest system 
leads to a significant complexation of the phase transi- 
tions. Phase transitions here are recognized by the lost 
of analyticity of the relevant thermodynamic potential V 
(Planck or Hclmholtz free energy, gran potential or any 
other admissible generalization). 

We have now five thermodynamical anomalies for the 
open systems in contract to the only two relevant when 
they are isolated. Besides the introduction of the zero- 
order phase transitions, we consider also a distinction 
among those continuous and discontinuous phase transi- 
tions which are related with a lost of analyticity of the 
entropy S. A tentative classification scheme is summa- 
rized in the Table ITT1 

Finally, the Table ITTT1 summarizes the "genealogy" of 
the phase transitions associated to the mathematical 
anomalies type A and B of the microcanonical entropy 
S which were described in the present section. While 



Thermodynamical anomaly 


Classification 


V is discontinuous as consequence 
of a discontinuity of control parame 
ter 77 = dS/dl used in the Legendre 
transformation V — nl — S . 


zero-order PT 


Some of first derivatives of V are 
discontinuous as consequence of 
the ensemble inequivalence. 


lst-kind 
discontinuous PT 


Some of first derivatives of V are 
discontinuous as consequence of 
a discontinuity in some of the first 
derivatives of the entropy S. 


2st-kind 
discontinuous PT 


Lost of analyticity of V with con- 
tinuous first derivatives without 
a lost of analyticity of entropy S. 


lst-kind 
continuous PT 


Lost of analyticity of V with con- 
tinuous first derivatives associated 
to a lost of analyticity of entropy S. 


2st-kind 
continuous PT 



TABLE II Tentative classification scheme of the phase tran- 
sitions (PT) within a "canonical description" (open system). 



Anomaly 


Isolated system 


Open system 


Type A 


microcanonically 
irrelevant 


lst-kind discont. PT 
lst-kind cont. PT 


Type B.I 


micro, continuous PT 


2st-kind cont. PT 


Type B.II 


micro, discontinuous PT 


zero-orden PT 
2st-kind discont. PT 



TABLE III Genealogy between the thermodynamic anoma- 
lies type A an B under different external conditions. 



the anomaly type A (non concavity of S in a given 
reparametrization 0) is irrelevant when the interest sys- 
tem is isolated, the incidence of certain experimental 
setup turns this behavior in a lost of analyticity of the rel- 
evant thermodynamical potential V associated with the 
ensemble inequivalence (a lst-kind discontinuous PT or 
its limiting case: lst-kind continuous phase transition). 
Type B anomalies always leads to the lost of analytic- 
ity of V whenever they are not hidden by the ensem- 
ble inequivalence. The corresponding phase transitions 
are always originated from the occurrence of ergodicity 
breaking. While microcanonical continuous phase transi- 
tions are directly related to the 2st-kind continuous phase 
transitions, a microcanonical discontinuous phase transi- 
tions can manifest as a zero-order phase transition or a 
2st-kind discontinuous phase transition. 
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VI. CONCLUSIONS 

We have shown in the section [H] that the microcanon- 
ical description is characterized by the presence of an in- 
ternal symmetry whose existence is closely related to the 
dynamical origin of this ensemble: the reparametrization 
invariance. Such symmetry leads naturally to a new geo- 
metric formalism of the Thermodynamics within the mi- 
crocanonical ensemble which is not based on the consid- 
eration of a Riemannian metric derived from the Hessian 
of the microcanonical entropy like other thermodynamic 
formalisms proposed in the past ijlffi . While the micro- 
canonical entropy becomes a scalar function within the 
present geometrical framework, we show that its convex 
or concave character depends on the reparametrization of 
the integrals of motion which are relevant for the micro- 
canonical description. Such ambiguity leads necessarily 
to a reconsideration of any classification scheme of the 
phase transitions based on the concavity of the micro- 
canonical entropy (@). 

Sections ITTtIvI were devoted to carry out a tentative 
analysis of the above question. Interestingly, this aim 
demands a necessary and unexpected generalization of 
the Gibbs canonical ensemble and the classical fluctua- 
tion theory where the reparametrization invariance in- 
troduced in the section ^ plays a more fundamental 
role than in the conventional Thermodynamics. While 
reparametrization changes do not alter the microcanon- 
ical description, such transformations represent specific 
substitutions of the external experimental setup leading 
to different generalized canonical descriptions 1281 for the 
open system. Particularly, we have shown that the well- 
known ensemble inequivalence of the Gibbs canonical en- 
semble^] can be avoided within an appropriate general- 
ized canonical ensemble 1281 which describes an open sys- 
tem put in contact with a generalized version of the Gibbs 
thermostat (heat bath) whose inverse temperature f3 ex- 
hibits correlated fluctuations with the fluctuations of the 
total energy E of the interest system (see in subsection 
IIILCf) . This feature of the generalized canonical ensem- 
ble put the basis for an unexpected generalization of the 
classical fluctuation theory (|13l) where the inverse tem- 
perature of the generalized thermostat [3 and the total 
energy of the interest system E behave as complementary 
thermodynamic quantities and the convex regions derived 
from the entropy Hessian admit a simple interpretation: 
there the inverse temperature (3 can not be fixed since 
(S(3SE) > 1 (see in subsection IIILE|) . The analysis of 
such questions demands a further study. Particularly, we 
have only consider in the present work a generalization 
of the Gibbs canonical ensemble, while it is reasonable a 
generalization of any Boltzmann-Gibbs distributions. 

The possibility of avoid the ensemble inequivalence by 
using the generalized canonical ensemble 1281 allows us to 
improve some Monte Carlo methods based on the consid- 
eration of the Gibbs canonical ensemble. This aim was 
carried out in the subsection IIII.DI and the ref. l(20|) as a 
example of application of the present reparametrization 



invariance ideas for the implementation of a generalized 
version of the well-known Metropolis importance sample 
algorithm ifT^ . 

The most important conclusion obtained from our sub- 
sequent analysis is the recognition about the fundamen- 
tal role of the arbitrary external conditions (the using of 
different experimental setups) which can be imposed to 
the interest system in order to control its thermodynamic 
equilibrium. Particularly, the existence and nature of the 
phase transitions (admitting them as a sudden change of 
the thermodynamic behavior of a given system while a 
smooth change of a control parameter) depends crucially 
on the nature of the external experimental setup. This 
fact was shown in the section Hvl bv considering the mi- 
crocanonical description of the q = 10 states Potts model 
system, where the Gibbs canonical description predicts 
the occurrence of a first-order phase transition and the 
microcanonical description exhibits two microcanonical 
phase transitions which are hidden by the ensemble in- 
equivalence. 

While there are only two thermodynamical anomalies 
which can be classified as phase transitions within the 
microcanonical ensemble (presumably associated to the 
occurrence of ergodicity breaking phenomenon (0)), there 
are five of them in an open system ("within a general- 
ized canonical description" ) : two of them closely related 
with the ensemble inequivalence (lst-kind discontinuous 
PT and its limiting case, the lst-kind continuous PT) 
and the other three with the lost of analyticity of the 
entropy S in the thermodynamic limit (zero-order PT, 
2st-kind discontinuous PT and the 2st-kind continuous 
PT). To our viewpoint the present classification scheme 
is only tentative and also deserves a further analysis of 
some paradigmatic models as well as the revision of the 
available literature. 



APPENDIX A: Riemannian metric from the Hes- 
sian of entropy 

Ruppcincr in ref. (ll3l) considers a subsystem a of a 
larger closed system £ in a thermodynamic equilibrium 
characterized by the macroscopic parameters xo ( a sys- 
tem + its environment - a situation schematically repre- 
sented in FIG llOfl . The probability that the subsystem a 
be in a the macrostate characterized by the parameters 
x, whose values are enclosed in the range (x, x + dx) is 
given by: 



P (x \x ) dx = —— — -Q,(x,x )dx, (Al) 
W(x ) 

being Q(x, Xq) the density of states at the configura- 
tion, and W (xq) = J 57 (x, xq) dx, the total number of 
microstates of the closed system. Instead of using the 
Boltzmann entropy associated to the closed system E, 
Sb(xo) = ]&W(xq), this author deals with the "en- 
tropy" So (x, x ) = In Q (x, x ) . 
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FIG. 10 Usual setup in classical thermodynamics fluctuation 
theory and the canonical ensembles in Statistical Mechanics: 
a subsystem a of a larger closed enviroment E, which are 
ordinarily regarded as uniform systems. 



S e = In {il (E, N) 5eq} where 5eo is a very small en- 
ergy constant which defines certain coarsed grained par- 
tition V e of the phase space T. The entropy in other 
reparametrization (E) = Nip (E/N) could be esti- 
mated by S v (6, N) = \n{fl(@,N)5ip } , where the 
small constant Sipo defines another coarsed grained par- 
tition Vip of the phase space. Obviously, the entropy 
estimates S e and S v are not equal since the partitions 
V € and V^p are essentially different. Taking into consid- 
eration the transformation rule 



de 

dE 



f2 (0, N) = 
the absolute discrepancy SS: 
SS = S V (Q,N)- S S (E,N) 



n(E,N), 



(Bl) 



depends on the relative energy e = E/N, but the relative 
discrepancy: 



The most probable macroscopic configuration x m of 
the subsystem a can be derived as a solution of the sta- 
tionary conditions: 



dS Q (x m ,x ) 



0, 



(A2) 



which ensures the nonnegative definition of the matrix 
g^y defined from the Hessian : 



6V ( x o) 



1 d 2 Sg (x m ,x ) 
V dx^dx" 



(A3) 



The matrix IA3I is invariant under " reparametrizations" 
of the thermodynamic parameters x of the subsystem a, 
y = y(x), since the second terms of the transformation 
rule: 

a 2 So (x (y) ,x ) _ dx» dx v d 2 S d 2 x^ dS 
dy a dy> 3 ~ dy a dyP dx»dx u dy a dyl 3 dx» ' 

dx^ 1 dx^ 

9ap (x ) = Q^Q^g^ (*o) (A4) 

vanish as a consequence of the stationary conditions IA2I 
This feature allows to take the matrix g^ v as a Rieman- 
nian metric which provides important information about 
the fluctuations around the most probable macroscopic 
configurations of the subsystem a. The reader can com- 
pare the result lA"4l with the one expressed in the equation 
1 151 in order to convince himself that the entropy Hessian 
of an isolated system is not a Riemannian tensor. 



APPENDIX B: Coarsed grained microcanonical 
entropy 

For practical purposes, the Boltzmann entropy Sb = 
In W can be estimated by the coarsed grained entropy 



ss_ 
~~s 



— with S = i [S e 



Srr 



(B3) 



drops to zero in the thermodynamic limit N — > oo as con- 
sequence of the extensive character of the entropy when 
N is large enough. 

The above reasonings allows us to consider the coarsed 
grained entropy, or more exactly, the entropy per particle 
in the thermodynamic limit: 



s (e) = lim 



S (eN, N) 
N 



= lim 

iV— >oc 



S (8 (eN) , N) 
N 



s(ip) 



as a continuous scalar function whose form is independent 
from the nature of the coarsed grained partitions V e and 
V v whenever their respective volume constants Seq and 
Sipo be small enough. Thus, we can assume along the 
present work the approximation S (E, N) = S (0, N). 



APPENDIX C: Others covariant quantities derived 
from the Boltzmann entropy 

Let us suppose that the integrals of motion depend on 
certain external parameter a, I a = I(X;a). By taking 
the partial derivative of the relation: 



n(0)= dXO(X)S I-I(X;a) 



we obtain the identity: 



f2 da 



{n<(M?>}, 



(Cl) 



(C2) 



where the microscopic quantity = F^ (X;a) repre- 
sents the k-th component of the generalized force associ- 
ated to the parameter a: 



81 s 
da 



(C3) 
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It is very easy to verify that the expectation values 

v s = (^OFf^ correspond to the components of a con- 

travariant vector because of they obey to the transfor- 
mation rule: 



JK 



d<p 



K 



di s 



The differential operator: 



V(u) 



1 d 



{nu s }, 



(C4) 



(C5) 



is just a linear transformation becoming a contravariant 
vector field u = {u s } in a scalar function. 

The invariant volume element dfj, g in the Riemannian 
geometry is derived from the metric g^ v as follows: dfi g = 
y/gd n x, where ^fg = -y/ |det <? M „|. The divergence of a 
given vectorial field u in this framework is given by: 



divu = {y/gu 1 *} . 



(C6) 



The analogy between the invariant volume elements 
d[i g = ■Jgd n x r~j dfj, = Qdl allows us to consider the 
differential operator T> (u) IC5l as the divergence of a con- 
travariant vector field. Thus, it is obvious that the ther- 
modynamic identity IC2I is reparametrization invariant. 

Substituting = 1 and O = in the identitv lU2l and 
denoting by S = In W where W = £151, being 51 a small 
constant volume, the following equations are straightfor- 
wardly obtained: 



dS_ 

da 



F 1 



dS_ 

da 
dF T 
da 



dS 
dS 



F„ 



F 



TS 



dFl 
dl s ' 
d 

f dF 7 



F 



TS 



(C7) 
(C8) 



where we use the following nomenclature for the micro- 
canonical averages F% = (.F^\ and Fj 5 = /fJF^V 

The second equation can be rephrased by introducing 
the correlation tensor G^ s = FT S — FjF„ s as follows: 



da 



dS 
dl s 



G, 



TS 



di s a a di s ' 



(C9) 



In the thermodynamic limit the quantities dFf /dl s and 
can be disregarded, and the quantity <S* be- 
comes in the microcanonical entropy obtained from the 
coarsed grained of the phase space. Thus, the thermody- 
namic identities: 



dS_ 

da 



d^ F s dFl 
di s a ' da 



di s a a di s: 



(CIO) 



express the relation among the entropy, the generalized 
forces, the response functions Xa = 9F^ /da and the cor- 
responding correlation tensor G„ associated to certain 



external parameter a within the microcanonical ensem- 
ble. Particular expressions of these results are the fol- 
lowing: 



dS_ 
dV 
dS_ 



0P, 



(3M 



dp 

dV = 
dM 

~dB 



/3G p +p 



dp 
dE' 



M 



dM 
~dE' 



(Cll) 
(C12) 



for a fluid (j> - pressure, V - volume, G p - pressure dis- 
persion) and a magnetic system (M - magnetization, B 
- magnetic field, Gm - magnetization dispersion) respec- 
tively, where (3 — dS/dE is the microcanonical inverse 
temperature. 
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